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Dynamic  Stability  Experiment  of  Maglev  Systems 

Y.  Cai,  D.  M.  Rote,  T.  M.  Mulcahy,  Z.  Wang,  S.  S.  Chen,  and  S.  Zhu 


Abstract 


This  report  summarizes  the  research  performed  on  maglev  vehicle  d5mamic 
stability  at  Argonne  National  Laboratory  during  the  past  few  years.  It  also 
documents  magnetic-force  data  obtained  from  both  measurements  and 
calculations.  Because  dynamic  instability  is  not  acceptable  for  any  commercial 
maglev  system,  it  is  important  to  consider  this  phenomenon  in  the  development  of 
all  maglev  systems.  This  report  presents  dynamic  stability  experiments  on 
maglev  systems  and  compares  their  numerical  simulation  with  predictions 
calculated  by  a  nonlinear  dynamic  computer  code.  Instabilities  of  an 
electrodynamic  system  (EDS)-type  vehicle  model  were  obtained  from  both 
experimental  observations  and  computer  simulations  for  a  five-degree-of-freedom 
maglev  vehicle  moving  on  a  guideway  consisting  of  double  L-shaped  aluminum 
segments  attached  to  a  rotating  wheel.  The  experimental  and  theoretical 
analyses  developed  in  this  study  identify  basic  stability  characteristics  and  future 
research  needs  of  maglev  systems. 


1  Introduction 


The  dynamic  response  of  maglev  systems  is  important  for  safety  and  ride 
quality,  guideway  design,  and  system  costs.  Ride  quality  is  determined  by  vehicle 
response  and  environmental  factors  such  as  humidity  and  noise.  The  dynamic 
response  of  vehicles  is  the  key  element  in  the  determination  of  ride  quality,  and 
vehicle  stability  is  one  of  the  important  elements  relative  to  safety.  To  design  a 
guideway  that  provides  acceptable  ride  quality  in  the  stable  region,  the  vehicle 
dynamics  must  be  understood.  Furthermore,  the  trade-off  between  guideway 
smoothness  and  levitation  and  control  systems  must  be  considered  if  maglev 
systems  are  to  be  economically  feasible.  Vehicle  dynamics,  which  link  the 
guideway  with  other  maglev  components,  must  be  analyzed  and  tested  in  detail 
when  considering  a  maglev  system  for  commercial  use  (Chen  et  al.,  1992;  Coffey 
et  al.,  1991). 
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For  safety,  maglev  systems  should  be  stable.  Because  dynamic  instabilities 
are  not  acceptable  for  any  commercial  maglev  system,  it  is  important  to  consider 
these  phenomena  when  designing  and  developing  maglev  suspension  systems. 
However,  to  date,  this  issue  has  not  received  much  attention  worldwide  (Masada, 
1993;  Wiescholek  et  al.,  1993;  Coffey,  1993;  Ozeki,  1993;  Kim  et  al.,  1993;  Martinelli 
and  Morini,  1993;  Lian  et  al.,  1993;  and  Wyczalek,  1993)  even  though  instabilities 
of  maglev  systems  were  confirmed  earlier  by  several  investigators  (Chu  and 
Moon,  1983;  Chen  et  al.,  1992;  Cai  et  al.,  1992b). 

Therefore,  the  objective  of  the  R&D  program  on  dynamic  stability  of  maglev 
systems  at  Argonne  National  Laboratory  (ANL)  during  the  past  four  years  (Chen 
et  al.,  1992;  Cai  et  al.,  1992a;  1992b;  1992c;  1993a;  1993b;  and  1993c;  Chen  et  al., 
1993;  Zhu  et  al.,  1994;  Cai  and  Chen,  1995)  has  been  to  develop  a  general  approach 
for  investigating  and  predicting  the  occurrence  of  instabilities  in  maglev 
suspension  systems  and  to  provide  a  better  understanding  of  the  conditions 
(design  features  and  parameter  values)  that  can  lead  to  d3mamic  instabilities. 

The  repulsive  levitation  system,  or  the  so-called  electrodynamic  system 
(EDS),  is  often  thought  to  be  inherently  stable.  However,  its  response  to 
perturbations  is  frequently  unstable  and  susceptible  to  catastrophic  oscillations, 
particularly  in  rectangular- trough  configurations.  So  far,  only  a  few  analytical 
and  experimental  studies  (Cai  et  al.,  1992a;  1992b;  1993a;  and  1993b)  have  been 
performed  to  gain  an  understanding  of  the  stability  characteristics  of  EDS-type 
maglev  systems.  Davis  and  Wilkie  (1971)  studied  a  magnetic  coil  moving  over  a 
conducting  track  and  concluded  that  negative  damping  occurs  at  velocities 
greater  than  the  characteristic  velocity  based  on  thin- track  theory.  Ohno  et  al. 
(1973)  studied  the  pulsating  lift  forces  in  a  linear  synchronous  motor.  These 
pulsating  forces  may  cause  parametric  and  combination  resonance,  in  addition  to 
heaving  and  pitching  oscillations.  Experiments  on  the  Massachusetts  Institute  of 
Technology  (MIT)  magneplane  showed  obvious  evidence  of  dynamic  instabilities 
on  film  in  the  early  1970s,  but  dynamic  stability  was  not  studied  in  detail.  An 
experimental  vehicle  with  three  degrees  of  freedom  (DOF),  floating  above  a  large 
rotating  wheel,  was  found  by  Moon  (1974)  to  have  a  lateral-roll-yaw  instability. 
Also,  experiments  performed  at  MIT  on  a  test  track  showed  pitch-heave 
instability.  Negative  magnetic  damping  was  demonstrated,  but  was  dominated  by 
aerodynamic  damping  (Moon,  1977).  A  conducting  guideway,  consisting  of  L- 
shaped  aluminum  segments  attached  to  a  rotating  wheel  to  simulate  the  full- 
scale  Japanese  guideway  at  Miyazaki,  was  studied  experimentally  and 
analytically  by  Chu  and  Moon  (1983).  Divergence  and  flutter  of  a  vehicle  model 
with  two  DOF  were  obtained  for  coupled  yaw-lateral  vibration;  the  divergence 
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leads  to  two  stable  equilibrium  yaw  positions,  and  the  flutter  instability  leads  to  a 
limit  cycle  of  coupled  yaw  and  lateral  motions  near  the  magnetic  drag  peak. 

Because  very  limited  studies  have  been  published  and  many  stability  issues 
remain  unresolved,  the  theory  of  motion-dependent  magnetic-force-induced 
instability  (which  consists  of  both  quasistatic  motion  and  unsteady  motion 
theories)  was  developed  in  our  previous  work  (Cai  et  al.,  1992a;  1992h;  1993a).  A 
method  of  obtaining  motion-dependent  magnetic-force  coefficients  experimentally 
was  identified  and  used  to  measure  quasistatic  motion-dependent  magnetic-force 
coefficients  of  a  maglev  system  with  a  double  L-shaped  guideway  from  an 
experimental  investigation  at  Argonne  National  Laboratory  (Cai  et  al.,  1992b). 
Motion-dependent  magnetic  forces  are  the  controlling  elements  in  the  stability  of 
maglev  systems.  Therefore,  this  earlier  integrated  analytical/experimental  study 
presented  a  systematic  method  to  determine  the  stability  of  maglev  systems  and 
answered  a  series  of  questions  on  maglev  stability.  With  a  simplified  analytical 
three-DOF  vehicle  model,  divergence  and  flutter  were  observed  for  coupled 
vibration  of  a  vehicle  on  a  guideway  consisting  of  double  L-shaped  aluminum 
segments  (Cai  et  al.,  1992a;  1992b).  Divergence  and  flutter  of  this  analytical  model 
were  later  verified  by  numerical  solution  of  a  case  with  the  same  vehicle 
configuration  (Cai  et  al.,  1993a;  1993b;  Cai  and  Chen,  1995). 

At  Argonne  National  Laboratory  (ANL),  an  analysis  of  dynamic  instabilities 
of  an  EDS-t3q)e  maglev  suspension  system,  with  three-  and  five-DOF  vehicles 
traveling  on  a  double-L  shaped  set  of  guideway  conductors,  was  conducted  (Cai  et 
al.,  1993a;  1993b;  Cai  and  Chen,  1995).  Both  analytical  and  numerical  approaches 
were  used,  and  various  magnetic  suspension  forces,  compiled  from  experimental 
data,  were  incorporated  into  the  theoretical  models.  Divergence  and  flutter  were 
obtained  from  analytical  and  numerical  solutions  for  coupled  vibration  of  the 
three-DOF  maglev  vehicle  model.  A  computer  code  for  numerically  simulating 
d3mamic  stability  of  the  five-DOF  vehicle  model  was  developed,  and  extensive 
computations  with  various  parameters  were  performed  to  determine  the  stability 
characteristics  of  EDS-type  maglev  systems.  Instabilities  of  five  directions  of 
motion  (heave,  slip,  roll,  pitch,  and  yaw)  of  the  dynamic  vehicle  model  were 
observed  and  it  was  demonstrated  that  system  parameters  such  as  system 
damping,  vehicle  geometry,  and  coupling  effects  among  five  different  motions 
play  very  important  roles  in  the  occurrence  of  dynamic  instabilities  in  maglev 
systems  (Cai  et  al.,  1993a). 

Magnetic  forces  are  basically  position-dependent,  although  some  are  also 
velocity-dependent.  These  motion-dependent  magnetic  forces  can  induce  various 
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types  of  instability.  In  addition,  the  periodic  structure  of  the  motion-dependent 
magnetic  forces  may,  in  some  cases,  also  induce  parametric  and  combination 
resonances  (Cai  et  al.,  1992b).  To  date,  magnetic  forces  are  obtained  mainly  from 
experimental  measurements.  At  ANL,  two  different  experiments  to  measure 
force  were  performed.  In  the  first  experiment,  a  quasistatic  method  was  applied 
to  measure  magnetic  forces  when  a  permanent  magnet  was  moving  over  an 
aluminum  L- shaped  ring  mounted  on  the  top  surface  of  a  1.2-m  diameter  rotating 
wheel.  The  measured  forces  are  a  function  of  position  at  a  given  steady  speed  (Cai 
et  al.,  1992h).  In  the  second  experiment,  d5mamic  (motion-dependent)  magnetic 
force  measurements  were  based  on  an  unsteady- motion  theory.  A  direct  method 
was  used  to  measure  magnetic  damping  and  stiffness,  and  to  investigate  the 
effect  of  various  parameters,  such  as  conductivity,  gap,  excitation  frequency,  and 
oscillation  amplitude  (Chen  et  al.,  1993;  Zhu  et  al.,  1994).  Experimental  results 
indicated  that  negative  magnetic  damping  will  develop  once  the  characteristic 
speed  is  exceeded.  In  the  maglev  system,  instability  will  occur  when  the  damping 
value  becomes  negative.  A  computer  model  simulation  of  magnetic  damping 
forces  in  a  maglev  system  confirmed  the  existence  of  negative  damping 
phenomenon  (He  et  al.,  1994). 

However,  in  this  report,  we  include  only  the  magnetic  force  data  measured  by 
the  quasistatic  method.  In  recent  years,  experiments  were  conducted  at  ANL  to 
measure  those  magnetic  forces.  We  summarize  those  magnetic  forces  in  this 
report  because  those  data,  obtained  from  steady  state  experiments,  play  a  very 
important  role  in  numerical  simulation  of  dynamic  stability  analysis  of  maglev 
vehicles.  All  of  the  experimental  data  have  heen  incorporated  into  an  existing 
computer  program  for  nonlinear- dynamics  simulation  of  maglev  vehicles.  We 
also  believe  that  integrating  magnetic  force  data  that  have  never  been  published 
will  provide  potential  references  for  future  maglev  programs. 

The  main  purpose  of  this  report  is  to  present  an  overview  of  recent  extensive 
experimental  investigation  on  the  dynamic  stability  of  maglev  systems  at  ANL. 
Two  series  of  tests  were  performed  with  a  free  vehicle  moving  on  a  double 
L-shaped  aluminum  guideway  mounted  on  the  top  of  a  rotating  wheel.  In  Test  1, 
a  vehicle  model  was  supported  hy  four  permanent  magnets  on  four  corners, 
whereas,  in  Test  2,  four  magnets  for  levitation  and  four  magnets  for  guidance 
were  attached  to  the  vehicle.  In  both  tests,  the  vehicle,  constrained  in  its 
longitudinal  direction  by  a  metal  tether  attached  to  its  front,  was  free  to  move  in 
up  to  five  modes  (vertical  heave,  lateral  slip,  pitch,  yaw,  and  roll).  Several 
accelerometers  and  a  force  transducer  (only  for  drag  direction)  were  placed  on  the 
vehicle.  The  vehicle's  motion  in  five  directions  was  measured  by  double  analog 
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integration.  When  the  rotating  speed  of  the  wheel  varied,  stable  and  unstable 
motions  of  the  maglev  vehicle  were  observed  and  recorded. 

A  computer  program,  incorporating  all  of  the  quasistatic  magnetic  forces 
that  we  measured,  was  developed  at  ANL  for  simulating  the  dynamic  response  of 
maglev  systems.  The  program  can  simulate  nonlinear  dynamic  response  of 
maglev  systems  with  six  DOF  when  the  user  inputs  vehicle  and  track 
configurations;  it  was  employed  in  this  study  to  calculate  the  dynamic  response 
and  onset  of  instability  in  a  maglev  vehicle  model  with  five  DOF.  The  simulations 
and  predictions  were  compared  with  results  obtained  from  dynamic  stability 
experiments. 


2  Measurement  of  Magnetic  Forces _ _ 

Magnetic  forces  are  needed  for  any  analysis  of  vehicle  dynamics,  design  of 
guideway  structures  and  fastening,  and  prediction  of  ride  quality.  These  force 
components  are  considered  from  the  standpoint  of  vehicle  stability. 

As  a  significant  part  of  ANL's  maglev  R&D  program,  magnetic  forces  of  a 
permanent  magnet  moving  over  an  aluminum  L-shaped  sheet  guideway  were 
measured  experimentally  with  a  quasistatic  method;  this  report  provides  an 
overview  of  that  effort.  Details  of  the  experimental  setup  and  measimed  force  data 
are  presented  in  this  section.  The  measured  magnetic  forces  are  compared  with 
magnetic  forces  calculated  by  various  commercial  computer  codes,  including 
MAGFORCE  and  ELEKTRA. 

We  integrated  the  experimentally  determined  magnetic  forces  in  this  report 
because  data  yielded  by  experiments  play  a  very  important  role  in  numerical 
simulation  of  dynamic  stability  analysis  of  maglev  vehicles.  All  of  the 
experimental  data  have  been  incorporated  into  an  existing  computer  program  for 
nonlinear-dynamics  simulation  of  maglev  vehicles. 


2.1  Experimental  Setup  for  Measuring  Magnetic  Forces 

An  experiment  was  conducted  at  ANL  to  investigate  the  lift,  drag,  and 
guidance  magnetic  forces  on  an  NdFeB  permanent  magnet  moving  over  an 
aluminum  (6061-T6)  L-shaped  ring  mounted  on  the  top  surface  of  a  1.2-m 
diameter  rotating  wheel  (shown  in  Figs.  1  and  2).  For  a  given  rotating  speed  of 


76.2  mm 


Fig.  2.  Experimental  apparatus  for  measuring  magnetic  force 


Fig.  3.  Schematic  diagram  of  apparatus  used  to  measure 
magnetic  forces  on  L-shaped  aluminum  sheet 
guideway 
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the  wheel  (ranging  from  0  to  600  rpm,  or  0  to  37.7  m/s),  the  lift,  guidance,  and 
drag  magnetic  forces  were  measured  while  the  guidance  gap  Y*  and  lift  height  h 
were  varied.  A  schematic  diagram  of  the  measurement  approach  is  shown  in 
Fig.  3  (Cai  et  al,  1992). 

During  testing,  the  long  side  of  the  25.4  x  50.8  x  6.35-mm  rectangular  magnet 
was  oriented  parallel  to  the  direction  of  motion  of  the  L-shaped  guideway  and  was 
held  stationary  by  a  two-component  force  transducer  that  comprised  two  BLH 
C2G1  load  cells  connected  in  series  to  measure  lift  and  guidance  forces 
simultaneously.  Laboratory  weights  were  used  to  calibrate  the  transducer  and  to 
assess  crosstalk  (which  was  found  to  be  <2%).  The  base  of  the  load  cell  assembly 
was  mounted  on  motorized  stages  that  provided  accurate  positioning  (+0.05  mm). 
Out-of-roundness  of  the  L-shaped  guideway  ring  varied,  but  was  always  <±0.15 
mm  for  the  lateral  leg  and  +0.35  mm  for  the  vertical  leg.  Ability  to  exactly  position 
the  magnet  with  respect  to  the  guideway  dominated  our  experimental  error, 
which  was  estimated  as  ±5%. 

Lift,  guidance  and  drag  forces  were  measured  as  a  function  of  h,  with  Y* 
fixed  at  5  and  12.7  mm,  or  as  a  function  of  Y*  with  h  fixed  at  7  and  12.7  mm  while 
the  surface  velocity  of  the  lateral  leg  of  the  guideway  varied  from  4.55  m/s  to 
36.1  m/s. 


2.2  Results  and  Discussion  of  Measured  Magnetic  Forces 

Measured  lift,  guidance,  and  drag  forces  are  plotted  in  Figs.  4,  5,  and  6, 
respectively,  as  a  function  of  h,  with  Y*  fixed,  or  as  a  function  of  Y*,  with  h  fixed, 
and  with  the  wheel  surface  velocity  =  4.55,  7.89,  10.62,  13.32,  18.20,  24.15,  30.10,  and 
36.10  m/s. 

The  following  characteristics  of  measured  lift,  guidance,  and  drag  forces 
were  noted  from  Figs.  4,  5,  and  6. 

•  The  qualitative  trends  in  the  lift,  guidance,  and  drag  force  data  for 
the  L-shaped  guideway  are  the  same  for  all  tested  velocities. 
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(b) 


(d) 


Fig.  4.  Measured  lift  magnetic  force,  for  various  wheel  surface  velocities  V,  as  a 
function  of  lift  height  when  (a)  Y*  =  5.0  mm  and  (b)  Y*  =  12.7  mm,  and  as 
a  function  of  gap  Y*  when  (c)  h  =  7.0  mm  and  (d)  h  =  12.7  mm 


With  fixed  gap,  the  lift  force  decreases  as  height  increases  (Figs.  4a 
and  4h);  with  fixed  height,  the  lift  force  remains  almost  the  same 
when  gap  varies  (Figs.  4c  and  4d).  With  fixed  gap,  the  guidance 
force  remains  the  same  (Figs.  5a  and  5h),  whereas  with  fixed 
height,  the  guidance  force  decreases  as  the  gap  increases  (Figs.  5c 
and  5d).  The  drag  force  decreases  as  either  gap  or  height  increases 
(Fig  6). 

Compared  with  the  lift  and  drag  forces,  the  value  of  the  guidance 
force  is  much  smaller  than  that  of  the  lift  and  drag  forces.  Without 
additional  guidance  magnets,  the  guidance  force  of  a  magnet  could 
not  maintain  the  vehicle  in  stable  motion. 


Guidance  Force,  N  Guidance  Force, 


10 


(b)  (d) 


Fig.  5.  Measured  guidance  magnetic  force,  for  various  wheel  surface  velocities 
V,  as  a  function  of  lift  height  when  (a)  Y*  =  5.0  mm  and  (b)  Y*  =  12.7  mm, 
and  as  a  function  of  gap  Y*  when  (c)  h  =  7.0  mm  and  (d)  h  =  12.7  mm 


•  Two  features  of  guidance  force  are  evident.  First,  a  maximum 
occurs  in  the  guidance  force  variation  with  respect  to  height 
variations  at  a  fixed  gap  (as  shown  in  Figs.  5a  and  5b),  a 
phenomenon  that  is  caused  by  the  corner  region  of  the  guideway. 
Second,  a  minimum  occurs  in  the  guidance  force  variation  with 
respect  to  gap  variations  at  a  fixed  height  (as  shown  in  Figs.  5c  and 
5d).  This  minimum  is  caused  by  the  edge  of  the  lateral  leg  of  the 
guideway.  As  shown  in  the  authors'  previous  publications,  the  first 
feature  is  associated  with  a  flutter  and  the  second  with  a  divergence 
instability  of  a  maglev  vehicle  model  with  three  DOF  (Cai  et  al., 
1992a  and  1992b). 
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(b)  (d) 

Fig.  6.  Measured  drag  magnetic  force,  for  various  wheel  surface  velocities  V,  as 
a  function  of  lift  height  when  (a)  Y*  =  5.0  mm  and  (b)  Y*  =  12.7  mm,  and 
as  a  function  of  gap  Y*  when  (c)  h  =  7.0  mm  and  (d)  h  =  12.7  mm 

The  measured  lift,  guidance,  and  drag  forces  are  replotted  in  Figs.  7,  8,  and 
9,  respectively,  as  a  function  of  velocity,  with  Y*  or  h  fixed.  Characteristics  of  the 
lift,  guidance,  and  drag  forces  from  Figs  7,  8,  and  9  can  be  summarized  as 
follows: 


When  velocity  increases,  the  lift  and  guidance  forces  overall 
increase  (Figs.  7  and  8);  whereas  drag  force  increases  in  the  lower 
velocity  range  and  decreases  in  the  higher  velocity  range  (Fig.  9). 

In  the  lower  velocity  range,  the  drag  force  is  quite  large  when 
compared  with  lift  and  guidance  forces;  therefore,  the  drag  force 
may  introduce  instability  into  a  maglev  vehicle  in  a  low  speed 
range. 
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Fig.  7.  Measured  lift  magnetic  force,  for  various  values  of  h  and  Y*,  as  a  function 
of  velocity  when  (a)  Y*  =  5.0  mm,  (b)  Y*  =  12.7  mm,  (c)  h  =  7.0  mm,  and 
(d)  h  =  12.7  mm 


The  magnetic  force  data  from  experimental  measurements  are  included  in 
the  Appendix  of  this  report.  The  curve  fits  to  magnetic  forces  in  Figs.  4-9  were 
derived  with  polynomial  expressions  and  input  into  a  computer  code  to  simulate 
the  dynamic  response  of  a  maglev  vehicle. 


2.3  Magnetic  Force  Computations 

A  series  of  computations  on  a  permanent  magnet  moving  over  an  L-shaped 
guideway  were  carried  out,  and  compared  with  experimental  magnetic  force 
data.  The  reason  for  conducting  force  computations  is  not  only  to  verify  experi¬ 
mental  force  data  but  also  to  provide  potential  theoretical  support  when  magnetic 
forces  cannot  be  obtained  from  experiments.  Therefore,  force  calculations  can 
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(b) 


(d) 


Fig.  8.  Measured  guidance  magnetic  force,  for  various  values  of  h  and  Y*,  as  a 
function  of  velocity  when  (a)  Y*  =  5.0  mm,  (b)  Y*  =  12.7  mm,  (c)  h  = 
7.0  mm,  and  (d)  h  =  12.7  mm 


be  employed  on  the  computer  program  of  the  maglev  system  to  analyze  vehicle 
dynamics.  The  force  computations  were  conducted  with  different  computer 
codes,  including  MAGFORCE  and  ELEKTRA,  with  the  configuration  and 
dimensions  of  the  magnet  and  L-shaped  guideway  as  shown  in  Fig.  3. 


2.3.1  MAGFORCE 

For  a  flat-plate  guideway  without  considering  transverse  edge  effects, 
MAGFORCE  works  very  well.  However,  for  an  L-shaped  guideway,  MAGFORCE 
produces  relatively  large  discrepancies  due  to  modeling  approximations;  a  finite 
thick  rectangular  permanent  magnet  is  represented  hy  a  single-turn  coil  and  the 
transverse  edge  of  the  guideway  is  ignored. 
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(a) 


(c) 


Fig.  9.  Measured  drag  magnetic  force,  for  various  values  of  h  and  Y*,  as  a 
function  of  velocity  when  (a)  Y*  =  5.0  mm,  (b)  Y*  =  12.7  mm,  (c)  h  =  7.0 
mm,  and  (d)  h  =  12.7  mm 


In  our  experiments,  we  suspended  a  permanent  magnet  over  an  L-shaped 
guideway.  Because  MAGFORCE  does  not  simulate  a  permanent  magnet  directly, 
an  equivalent  coil  must  be  used  to  simulate  the  permanent  magnet.  We  assumed 
that  the  dimensions  of  the  coil  are  the  same  as  those  of  the  magnet  (5.08  x 
2.54  cm),  and  located  the  coil  in  the  middle  of  the  magnet,  as  shown  in  Fig.  10. 
This  kind  of  coil  model  was  used  in  flat-plate  guideway  problems,  and  it  worked 
very  well.  However,  for  L-shaped  guideway  problems,  large  errors  were  produced, 
especially  for  drag  and  guidance  forces.  Figures  11-15  show  the  force  results, 
which  show  that  lift  forces  match  experimental  data  reasonably  well.  But 
discrepancies  up  to  50%  were  noted  for  drag  and  guidance  forces. 
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Fig.  10.  Magnet  and  equivalent  coil 
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Fig.  11.  Measured  and  MAGFORCE -calculated 
lift  forces  when  Y*  =  12.7  mm 


Fig.  12.  Measured  and  MAGFORCE -calculated 
guidance  forces  when  Y*  =  12.7  mm 


Fig.  13.  Measured  and  MAGFORCE -calculated 
drag  forces  when  Y*  ■=  12.7  mm 


Fig.  14.  Measured  and  MAGFORCE -calculated 
lift  forces  when  Y*  -  5  mm 
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Fig.  15.  Measured  and  MAGFORCE -calculated 
drag  forces  when  Y*  =  5  mm 
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A  possible  reason  for  the  discrepancies  is  that  the  coil  was  located  in  the 
middle  plane  of  the  magnet  instead  of  on  the  bottom  surface  of  the  magnet.  An 
alternative  was  proposed  that  it  should  also  be  located  elsewhere  but  not  on  the 
surface  facing  the  guidance  plate.  Because  the  thickness  of  the  magnet  is 
6.35  mm,  the  dimensions  of  equivalent  coils  are  assumed  to  be  4.45  x  1.91  cm,  as 
shown  in  Fig.  16.  MAGFORCE  results  obtained  with  this  alternative  coil  are 
shown  in  Figs.  17-19.  Discrepancies  are  smaller  than  those  of  the  previous  model, 
but  the  shapes  of  the  guidance  force  curves  obtained  from  MAGFORCE  are 
different  from  those  that  were  obtained  from  experiments. 

MAGFORCE  was  developed  at  ANL  by  H.  Coffey  and  is  based  on  a  published 
Fourier  transformation  formulation  (Coffey  et  al.,  1972). 

For  design  purposes,  MAGFORCE  is  useful  because  it  will  give  us  an 
approximation  of  the  steady-state  capabilities  for  levitation,  guidance,  and 
propulsion.  However,  it  does  not  model  a  permanent  magnet  and  L-shaped 
guideway  accurately  enough  to  predict  dynamic  effects,  which  also  depend  on 
force  gradients.  However,  it  is  not  accurate  enough  for  theoretical  analysis. 


2.3.2  ELEKTRA 

The  velocity  version  of  ELEKTRA  is  a  finite-element  computer  code, 
developed  by  Vector  Fields,  Ltd.;  it  is  capable  of  solving  three-dimensional  eddy- 
current  problems  with  a  moving  conductor. 

Theoretically,  ELEKTRA  can  be  used  to  analyze  any  kind  of  eddy  current 
problems  with  complex  geometries.  It  has  been  successfully  used  to  analyze  flat 
plate  guideway  and  edge  effect  problems.  Here,  we  use  it  to  analyze  the  L-shaped 
guideway  problem. 
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Fig.  16.  Alternative  magnet  and  equivalent  coil 
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Fig.  17.  Measured  and  MAGFORCE -calculated 
lift  forces  when  Y*  =  12. 7  mm 


Fig.  18.  Measured  and  MAGFORCE -calculated 
guidance  forces  when  Y*  =  12. 7  mm 
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Fig.  19.  Measured  and  MAGFORCE -calculated 
drag  forces  when  Y*  =  12. 7  mm 
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Since  ELEKTRA  does  not  simulate  a  permanent  magnet  directly,  an 
equivalent  coil  must  be  used  to  simulate  the  permanent  magnet.  Therefore, 
equivalent  ampere  turns  (A-T)  of  the  coil  must  be  determined.  One  way  to 
determine  the  ampere  turns  is  to  compare  the  amplitude  and  shape  of  the 
magnetic  field  in  the  gap  between  the  magnet  (coil)  and  the  conductor  sheet. 
When  this  technique  is  used,  a  coil  equivalent  to  the  magnet  used  in  the 
experiments  would  have  5,994  A-T.  The  outer  dimensions  of  the  coil  used  in 
ELEKTRA  correspond  exactly  to  the  outer  dimensions  of  the  permanent  magnet 
used  in  the  experiment,  and  the  coil  cross  section  is  1  mm  wide  by  6.35  mm  high. 

During  construction  of  the  mesh,  one  should  make  the  mesh  distribution 
around  and  inside  the  coil  as  fine  as  possible  to  obtain  more  accurate  force 
computations.  The  total  number  of  nodes  for  our  case  is  -13,000-16,000,  and 
computational  time  for  one  velocity  is  -2-8  h,  depending  upon  the  velocity. 

Figure  20  shows  lift,  drag,  and  guidance  forces  vs.  velocity  calculated  for  h  = 
12.7  mm  and  L  =  12.7  mm.  ELEKTRA  results  agree  well  with  experimental  data. 
However,  ELEKTRA  is  too  time  consuming.  It  sometimes  required  6-8  h  (2000- 
3000  iterations)  at  a  SUN  work  station  to  obtain  converged  solutions. 

Figures  21-23  illustrate  forces  vs.  velocity  for  L  =  5  mm.  Figures  24-26  show 
forces  vs.  velocity  for  h  =  7  mm.  All  ELEKTRA  results  are  in  good  agreement  with 
experimental  data.  However,  ELEKTRA  has  two  shortcomings:  it  is  too  time 
consuming,  as  mentioned  above,  and  it  sometimes  has  difficulty  in  achieving 
convergence,  even  though  the  velocity  is  not  high,  e.g.,  25  m/s.  In  addition. 


Fig.  20.  Measured  and  ELEKTRA-calculated  lift, 
guidance,  and  drag  forces  when  Y*  = 
12.7  mm 
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Fig.  21.  Measured  and  ELEKTRA-calculated 
lift  forces  when  Y*  =  5  mm 


Fig.  22.  Measured  and  ELEKTRA-calculated 
guidance  forces  when  Y*  =  5  mm 


Fig.  23.  Measured  and  ELEKTRA-calculated 
drag  forces  when  Y*  =  5  mm 
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Fig.  24.  Measured  and  ELEKTRA-calculated 
lift  forces  when  h  =  7  mm 


Fig.  25.  Measured  and  ELEKTRA-calculated 
guidance  forces  when  h  =  7  mm 


Fig.  26.  Measured  and  ELEKTRA-calculated 
drag  forces  when  h  =  7  mm 
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ELEKTRA  has  been  verified  as  a  means  of  calculating  steady-state  magnetic 
forces,  but  dynamic  forces  (of  damping)  must  still  be  measured. 


2.4  Measuring  Magnetic  Forces  of  Guidance  Magnets 

Because  the  25.4  x  50.8  x  6.35-mm  rectangular  magnet  shown  in  Fig.  3 
cannot  provide  sufficient  guidance  force  for  a  maglev  vehicle  moving  on  an 
L-shaped  aluminum  guideway,  additional  guidance  magnets  are  needed  to 
increase  guidance  force  and  vehicle  stability.  A  rectangular  12.7  x  50.8  x  6.35-mm 
magnet  was  selected  as  a  guidance  magnet.  The  measurement  of  lift  and  drag 
forces  was  conducted  by  setting  the  magnet  on  a  rotating  drum  covered  with  a 
14.61-cm-wide,  6.35-mm-thick  aluminum  sheet  rolled  to  a  diameter  of  26.99  cm. 
The  rotating  speed  can  be  varied  from  0  to  3500  rpm;  therefore,  the  speed  of  the 
wheel  is  0-50  m/s  (the  detailed  experimental  setup  can  be  found  in  Chen  et  al., 
1993).  The  configuration  of  the  guidance  magnet  is  shown  in  Fig.  27,  where  the 
distance  between  the  magnet  and  the  aluminum  sheet  is  h. 

The  measured  lift  and  drag  forces  of  the  guidance  magnet  as  a  function  of 
velocity  are  plotted  in  Fig.  28,  with  h  fixed  at  8.5,  13.5,  and  18.5  mm.  Lift  and  drag 
force  data  are  included  in  the  Appendix.  The  curve  fits  to  magnetic  forces  in 
Fig.  28  were  derived  with  polynomial  expressions  and  incorporated  into  a 
computer  code  for  dynamic  simulation  of  a  maglev  vehicle. 


magnet 


Fig.  27.  Guidance  magnet  on  aluminum  sheet  guideway 


23 


Fig.  28.  Measured  lift  and  drag  magnetic 
forces  of  guidance  magnet  as  a 
function  of  velocity 


3  Dynamic  Stability  Experiments _  ■ _ 

The  experimental  apparatus  for  maglev  dynamic  stability  investigation  is 
shown  in  Fig.  29.  Two  L-shaped  aluminum  tracks,  like  an  inverted  T,  are 
mounted  over  a  rotating  wheel.  The  maglev  vehicle  is  magnetically  levitated  over 
the  aluminum  track  when  the  wheel  rotates  with  constant  but  adjustable  speed. 
The  vehicle,  constrained  in  its  longitudinal  direction  by  a  metal  tether  attached  to 
its  front,  is  free  to  move  in  up  to  five  modes  (vertical  heave,  lateral  slip,  pitch,  yaw, 
and  roll).  The  wheel  diameter  is  1.2  m.  The  wheel  rotating  speed  can  be  adjusted 
from  0  to  600  rpm  (or  0  to  37.7  m/s  on  the  surface  of  the  wheel). 

A  force  transducer  is  attached  between  the  tether  and  support  frame; 
therefore,  the  drag  force  applied  to  the  vehicle  can  be  measured.  Several  small 
accelerometers  are  installed  on  the  vehicle  body  to  measure  vehicle  motion  by 
double  analog  integration.  Displacement  and  force  signals  are  first  filtered  by 
band-pass  filters  to  eliminate  low-  and  high-frequency  noises  and  then  digitized 
and  stored  in  the  analyzer.  These  signals  are  then  analyzed  to  obtain  frequencies 
and  displacements  of  vehicle  motion  as  a  function  of  wheel  speed. 


Fig.  29.  Schematic  diagram  of  experimental  apparatus  for 
testing  maglev  dynamic  stability 

For  safety,  the  experimental  apparatus  is  placed  in  an  underground  room 
with  heavy  concrete  covers.  A  video  camera  is  set  beside  the  wheel.  The  test  runs 
can  be  observed  from  a  TV  set  and  recorded  in  a  VCR  aboveground. 

A  fishing  line  is  attached  to  each  of  the  vehicle's  four  corners  to  control  the 
position  of  the  vehicle,  like  moving  a  puppet.  When  the  vehicle  is  levitated,  the 
fishing  lines  are  loose;  if  instability  occurs,  the  vehicle  must  be  pulled  out  of  the 
track  to  prevent  damage  to  the  vehicle  and  magnets;  when  the  experiments  are 
over,  the  vehicle  is  hung  over  the  wheel. 

The  object  of  this  experiment  was  to  investigate  the  existence  of  instability  in 
maglev  systems.  The  term  instability  refers  to  motions  that  depart  from 
equilibrium  state  and  are  distinguished  from  random  oscillations  due  to  track 
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irregularities  and  from  resonance-driven  oscillations  associated  with  the  circular 
nature  of  the  rotating  guideway  (Chu  and  Moon,  1983). 

Two  series  of  tests  (A  and  B)  were  performed  with  a  free  vehicle  moving  on  a 
double  L-shaped  aluminum  guideway  mounted  on  the  top  of  the  rotating  wheel. 
In  Test  A,  a  vehicle  model  was  supported  by  four  permanent  magnets  on  four 
corners,  whereas  in  Test  B,  four  magnets  for  levitation  and  four  magnets  for 
guidance  were  attached  to  the  vehicle.  The  vehicle's  motion  in  five  directions  was 
measured  by  double  analog  integration.  When  the  rotating  speed  of  the  wheel  was 
varied,  stable  and  unstable  motions  of  the  maglev  vehicle  were  observed  and 
recorded. 


3.1  Test  A 

A  vehicle  model  with  four  25.4  x  50.8  x  6.35-mm  rectangular  levitation 
magnets  is  shown  in  Fig.  30.  The  vehicle  body  was  made  of  fiberglass/epoxy  (GIO) 
sheet.  Dimensions  of  the  vehicle  are  indicated  in  Fig.  30.  A  shaft,  which  is 
attached  to  the  tether  and  placed  in  the  lateral  center,  can  move  back  and  forth  in 
the  longitudinal  direction,  which  constrains  the  DOF  of  the  vehicle.  Two 
accelerometers  are  installed  on  the  vehicle  to  measure  vertical  and  lateral 
motion. 

To  identify  the  regions  of  stability  of  the  vehicle  model,  several  runs  were 
performed  with  various  vehicle  conditions,  such  as,  changing  lateral  gaps  to  vary 
guidance  forces,  adding  mass  to  vary  vertical  gaps  and  lift  forces,  and  centering 
or  offsetting  the  shaft.  Dynamic  instability  was  observed  in  most  of  vehicle 
conditions. 

Figure  31  shows  typical  experimental  results.  In  this  test,  the  total  mass  of 
the  vehicle  was  2.12  kg.  The  lateral  gap  was  set  at  155.6  mm  (see  Fig.  30).  The  test 
was  run  from  a  high  speed  of  595  rpm  to  a  low  speed  of  128  rpm  (if  the  speed  is  too 
low,  the  lift  forces  cannot  levitate  the  vehicle).  Vertical  and  lateral  displacements 
were  measured  with  accelerometers  by  analog  integration.  Drag  force  was 
measured  by  a  force  transducer  between  the  tether  and  the  frame.  Three  cases, 
with  various  shaft  positions  (in  the  center,  and  offset  50.8  and  76.2  mm)  were 
plotted  in  Fig.  31,  in  terms  of  RMS  vertical  and  lateral  displacements  and  RMS 
force  as  function  of  wheel  speed.  At  least  four  speed  regions  were  found  by  both 
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A 


Fig.  30.  Vehicle  model  with  four  levitation  magnets  in  Test  A 


mm 
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measurement  and  observation.  For  example,  consider  the  case  with  the  shaft  in 
the  center  (Case  1).  (1)  When  the  speed  is  lower  than  132  rpm,  the  vehicle  is 
unstable,  with  high  frequency  flutter  due  to  a  large  drag  force.  (2)  When  the  speed 
is  between  132  and  245  rpm,  the  vehicle  is  stable,  amplitudes  of  displacement  and 
force  are  relatively  small,  and  the  power  spectral  densities  (PSD)  show  clear 
harmonic  peaks.  (3)  When  the  speed  is  between  245  and  471  rpm,  the  vehicle  is 
unstable  and  dominated  by  slide  and  yaw  instability  (divergence),  and  the 
oscillating  frequency  is  quite  low.  (4)  When  the  speed  is  over  47 1  rpm,  the  vehicle 
is  unstable  in  almost  every  direction  and  the  oscillation  amplitudes  in  five 
directions  increase  significantly.  Cases  2  and  3  show  the  same  trends  for  four 
speed  regions. 

From  observation  and  analysis,  the  vehicle  model  with  four  magnets  (see 
Fig.  30)  is  likely  to  develop  a  divergence  instability  because  the  guidance  force 
provided  by  the  magnet  on  the  L-shaped  guideway  is  quite  small  when  compared 
with  lift  and  drag  forces  (see  Section  2).  Therefore,  a  new  vehicle  with  additional 
guidance  magnets  was  assembled  and  tested  in  Test  B. 


3.2  Test  B 

A  vehicle  model  with  four  25.4  x  50.8  x  6.35-mm  levitation  magnets  and  four 
12.7  X  50.8  X  6.35-mm  guidance  magnets  is  schematically  shown  in  Fig.  32.  The 
clearances  between  the  sheet  guideway  and  guidance  and  levitation  magnets  can 
be  set  at  several  values.  Four  small  wheels  are  attached  to  the  vehicle  to  prevent 
damage  from  dynamic  instability.  The  vehicle  weight  is  2.502  kg  or  24.52  N. 
Moments  of  inertia  from  measurement  are  0.4389,  0.2840,  and  0.1920  kgm^  for  x, 
y,  and  z  axes,  respectively. 

In  addition  to  the  force  transducer,  six  accelerometers  were  placed  on  the 
vehicle  (Fig.  33);  four  (Zl,  Z2,  Z3,  and  Z4)  at  the  four  comers  of  the  vehicle  in  the 
vertical  direction  and  two  (Y 1  and  Y2)  at  two  comers  of  one  side  of  the  vehicle,  in 
the  lateral  direction.  Accelerations  were  doubly  integrated  to  obtain 
displacement.  Therefore,  vehicle  motion,  including  vertical  heave,  lateral  slide, 
pitch,  yaw,  and  roll  can  be  calculated  from  measured  displacements  as  follows: 


101.6  mm 
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Fig.  33.  Accelerometer  arrangement  for 
vehicle  model  in  Test  B 


Heave  (t)  =  [Zi  (t)  +  Z2  (t)  +  Z3  (t)  +  Z4  (t)]  /  4 

Slide  (t)  =  [Yi(t)  +  Y2(t)]/2 

Pitch  (t)  =  [Z3  (t)  +  Z4  (t)  -  Zi  (t)  -  Z2  (t)]  /  (2Li ) 

Yaw  (t)  =  [Y2  (t)  -  Yi  (t)]  /  L2 

RoU  (t)  =  [Zi  (t)  +  Z4  (t)  -  Z2  (t)  -  Z3  (t)]  /  (2W) 

Vehicle  response  and  stability  were  tested  as  a  function  of  speed  for  a  specific 
configuration  of  levitation  and  guidance  magnets.  When  the  flywheel  was 
running  from  300  to  600  rpm,  two  dynamic  instability  regions  were  noted.  Two 
video  tapes  showing  the  dynamic  response  of  the  vehicle  at  different  speeds  are 
available. 

A  test  was  conducted  to  see  if  the  vehicle  can  be  operated  in  a  stable  position 
without  the  four  guidance  magnets.  The  result  showed  that,  in  this  setup,  the 
vehicle  was  always  unstable. 
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Several  series  of  tests  were  performed  to  investigate  the  effect  of  the  clearance 
between  the  sheet  guideway  and  the  magnets.  The  gap  between  the  sheet 
guideway  and  the  levitation  magnets  was  set  at  1.27  and  2.54  mm,  while  the  gap 
between  the  sheet  guideway  and  the  guidance  magnets  was  set  at  8.5,  13.5, 
18.5  mm.  The  response  characteristics  depend  on  the  gaps  between  the  guideway 
and  levitation  or  guidance  magnets.  However,  some  general  characteristics  were 
noted. 

•  For  speeds  over  450  rpm,  the  motions  of  the  vehicle  were  small. 

•  For  speeds  between  350  and  450  rpm,  the  motions  of  the  vehicle  were 
fairly  large,  with  one  of  its  wheels  frequently  touching  the 
guideway. 

•  For  speeds  between  300  and  350  rpm,  the  motions  of  the  vehicle  were 
again  relatively  small. 

•  For  speeds  lower  than  300  rpm,  the  vehicle  was  not  yet  completely 
levitated. 

Figure  34  shows  results  of  a  detailed  test  performed  for  the  following  gaps: 
2.54  cm  between  the  vertical  part  of  the  guideway  L  and  the  edge  of  the  levitation 
magnets;  1.35  cm  between  the  vertical  part  of  the  guideway  L  and  guidance 
magnets  (see  Fig.  32).  Figure  34a  shows  RMS  longitudinal  force  as  a  function  of 
the  rotating  speed  of  the  wheel;  Fig.  34b  shows  RMS  displacements  at  six 
measurement  positions  as  a  function  of  speed;  pitch,  roll,  yaw,  heave,  and  slide 
motions  calculated  from  displacements  are  given  in  Fig.  34c  and  34d.  The  motion 
of  the  vehicle  at  various  rotating  speeds  was  recorded  on  a  video  tape  from  =600  to 
290  rpm.  At  each  of  the  following  speeds,  the  motions  were  recorded  for  about 
30  s:  596, 455,  440,  390,  380, 370, 330,  320,  310,  and  290  rpm. 

From  Fig.  34  and  recorded  video  tape  playback,  dynamic  instability  for  this 
vehicle  model  indeed  exists  and  large  motion  occurs  at  between  350  and  450  rpm. 
From  the  video,  it  was  noted  that  the  vehicle  touched  the  guideway  intermittently 
from  380  to  440  rpm;  while  at  370  and  455  rpm,  no  impact  was  noted.  The 
response  characteristics  depend  on  the  rotating  speed  of  the  flywheel.  At  440  rpm, 
the  motion  was  fairly  steady  with  regular  and  sometimes  intermittent  impacts. 
At  390  rpm,  the  motion  was  more  irregular;  it  appeared  to  be  chaotic  vibration.  At 
380  rpm,  the  motion  was  fairly  steady  but,  in  each  cycle,  there  were  several 
impacts. 
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Fig.  34.  Experimental  results  from  Test  B.  (a)  RMS  force;  (b)  RMS  displacement; 

(c)  degree  of  pitch,  roll,  and  yaw;  and  (d)  heave  and  lateral 
displacement,  all  as  a  function  of  velocity 


Significantly,  it  was  noted  that  slide  motion  was  much  smaller  in  Test  B  than 
in  Test  A  because  the  guidance  magnets  in  Test  B  provided  sufficient  guidance 
force.  From  experimental  data  and  observation,  oscillations  of  the  vehicle  were 
dominated  by  heave  and  roll  motions. 


Based  on  the  experimental  data,  it  is  difficult  to  identify  precisely  the  cause  of 
large-amplitude  oscillations.  Several  mechanisms  may  be  involved:  periodic 
excitation  because  of  imperfections  in  the  flywheel;  motion-dependent  magnetic 
forces;  aerodynamic  forces.  The  results  from  the  experiment  illustrate  future 
research  needs  in  the  area  of  vehicle  dynamics  of  maglev  systems.  To  predict 
vehicle  response  confidently,  magnetic  and  aerod3rnamic  forces,  and  guideway 
characteristics  must  be  quantified.  The  most  critical  elements  in  predicting 
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vehicle  response  are  magnetic  forces;  these  include  steady  and  unsteady 
magnetic  forces.  At  this  time,  only  limited  steady  magnetic  forces  are  available. 
Future  research  should  be  extended  to  the  prediction  and  measurement  of 
unsteady  magnetic  forces. 


4  Numerical  Simulation  of  Maglev  Vehicles  with  Five  DOF 

A  computer  code  for  simulating  the  dynamic  response  of  maglev  systems 
was  developed.  The  code,  which  incorporated  all  of  the  quasistatic  magnetic 
forces  that  were  measured  experimentally,  can  simulate  the  nonlinear  dynamic 
response  of  maglev  systems  with  six  DOF  when  the  user  inputs  vehicle  and  track 
configurations.  It  was  employed  in  this  study  to  calculate  dynamic  response  and 
the  onset  of  instability  in  a  maglev  vehicle  model  with  five  DOF.  The  simulations 
and  predictions  were  compared  with  results  obtained  from  dynamic  stability 
experiments. 


4.1  Theoretical  Aspects 

Consider  a  maglev  vehicle  system  with  six  DOF,  i.e.,  three  translations,  Ux, 
Uy,  and  Uz,  and  three  rotations,  0)x,  tOy,  and  (Oz,  as  illustrated  in  Fig.  35. 

Governing  equations  of  motion  of  this  six-DOF  maglev  system  are  derived 
from  Newton's  and  Euler's  equations,  and  the  transform  between  the  vehicle  and 
guideway  coordinate  systems  is  based  on  Euler  angles  (Cai  and  Chen,  1993a  and 
1995;  Coffey  et  al.,  1991). 

Consider  a  rigid  body  vehicle  with  six  DOF  (three  translations,  Ux,  Uy,  and 
Uz,  and  three  rotations,  CDx,  cOy,  and  coz-  Let  U  be  the  translation  vector  and  O  the 
rotation  vector,  i.e.. 


u  =  - 

Ux 

% 

0)x 

COy  > 

.Uz. 

COz 

(1) 
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Fig.  35.  Displacement  components  of  a 
maglev  system 


To  transform  from  the  vehicle  to  the  guideway  coordinate  system  (Coffey  et 
al.,  1991),  the  Euler  angles  y  (heading),  0  (attack),  and  0  (roll)  are  introduced.  A 
vector  Rg  in  the  guideway  coordinates  transforms  from  a  vector  Ry  in  the  vehicle 
coordinates  as 

Rg  =  MRv,  (2) 

with  the  transformation  matrix  given  by 

'cos6cosv|f  -sin\}/cos<})  +  cosy  sin9sin<|)  sinv|/sin(t)  +  cos\}rsin6cos({> 

[/]  =  cos0sinv|/  cosv|/cos0  +  sin\|/sin0sin4)  -cos\(/sin4)  +  sin\|/sin0cos^  .  (3) 
— sin0  cos  0  sin  4)  cos  0  cos  (j) 

Equations  of  motion  expressed  in  the  vehicle  coordinate  system  can  be 
described  as 
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MU  +  CU  =  MG  +  m(u  X 
lC  +  DtJ  =  ir  +  RFm 


(4) 


where  M  and  I  are  the  mass  of  the  vehicle  and  the  moment  of  inertia  about  the 
center  of  the  mass  inertia  moment  of  the  vehicle,  and  are  defined  hy 
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and  C  and  D  are  damping  coefficient  matrices 
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The  damping  coefficient  in  Eq.  6  should  be  determined  from  experimental  data. 
Before  these  experimental  data  are  available,  assumed  damping  ratios  can  be 
utilized  in  the  simulations.  The  elements  of  G  and  F  are  defined  by 
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Fjjj  is  the  magnetic  force  applied  on  the  center  of  the  mass,  and  R  is  the  moment 
arm  matrix.  Because  the  applied  forces  depend  on  the  distances  between  the 
vehicle  and  the  guideway  surfaces,  but  all  components  of  the  forces  and  moments 
in  the  equations  of  motion  are  to  be  expressed  along  the  axes  of  the  vehicle,  the 
transformation  must  convert  forces  from  one  coordinate  to  the  other. 


Because  the  time  derivatives  of  the  Euler  angles  \|/,  0,  and  <|)  are  related  to  the 
rotational  velocities. 


V  =  sec  0(0)2  cos  (|)  +  cby  sin  (!>) 

0  =  0)y  cos(l)-(b2sin(t) 

^  =  0)x  +  tan  0(0)2  cos  (])  +  d)y  sin  <|)) 


(9) 


Because  only  a  five-DOF  vehicle  model  is  considered  in  the  djmamic 
simulation  and  experiments  with  the  ANL  rotating-wheel  facility,  the 
longitudinal  motion  of  the  vehicle  Ux  is  restricted  and  can  be  assumed  to  be  zero  in 
the  equations  of  motion. 


4.2  Computer  Code  for  Maglev  Systems 

A  computer  code,  based  on  governing  Eqs.  4  to  9,  was  developed  for  nonlinear 
d5mamic  simulation  of  maglev  systems.  The  program  uses  the  fourth-order 
Runge-Kutta  method  to  numerically  integrate  nonlinear  differential  equations 
that  characterize  the  dynamics  of  maglev  vehicles. 

The  objective  of  the  computer  code  is  to  predict  the  dynamic  behavior  of  the 
five-DOF  maglev  systems  described  in  Section  3,  even  though  the  code  can  be 
extended  to  cover  a  wider  range  of  vehicle  and  guideway  designs.  The  different 
system  descriptions  in  Test  A  and  Test  B,  such  as  design  of  vehicle  and  track 
geometry,  number  and  size  of  magnets,  and  vehicle  traveling  speed,  can  easily  be 
included  in  a  user-friendly  input  package  of  the  code.  Then,  a  specific  subroutine 
of  the  code  will  automatically  calculate  essential  parameters  of  vehicle/track 
systems  for  later  numerical  simulation. 
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Magnetic-force  data  from  experimental  measurements  are  incorporated  in 
force  subroutines  of  the  code.  Once  the  distances  between  the  magnets  and  the 
sheet  guideway  are  calculated  out  during  each  time  step,  lift,  guidance,  and  drag 
forces  on  each  magnet  are  determined  from  prestored  magnetic-force  data  curves 
by  a  multidimensional  interpolation  method.  Then,  total  magnetic  forces  and 
moments  acting  on  the  vehicle  body  are  calculated  and  inserted  into  differential 
equations.  The  coordinate  transformation  (Eq.  2)  is  required  to  convert  forces 
from  one  coordinate  to  the  other,  because  the  applied  forces  are  dependent  on  the 
distance  between  the  vehicle  and  the  guideway  surfaces,  but  all  components  of  the 
forces  and  moments  in  the  equations  of  motion  are  to  be  expressed  along  the  axes 
of  the  vehicle  body. 

One  difficulty  with  describing  the  system  is  how  to  determine  system 
damping  in  the  differential  equation.  The  system  damping  is  very  important  in 
the  dynamic  and  stability  analysis.  However,  so  far,  an  actual  system  damping 
ratio  cannot  be  obtained  by  either  experimental  or  analytical  methods.  Therefore, 
we  assumed  the  reasonable  value  of  +2.5%  for  the  damping  ratio  in  most  of  the 
simulations.  Here,  the  reasonable  value  means  that,  with  this  value,  simulation 
results  would  be  closest  to  those  obtained  from  experiments.  Also,  different 
damping  values  are  considered  in  the  simulation.  Usually,  if  the  damping  ratio 
is  increasing  to  a  certain  value,  the  vehicle  motion  would  be  stable  in  the  whole 
speed  range;  if  less  damping  is  assumed,  the  instability  region  broadens.  The 
system  damping  affects  motions  in  all  directions. 


4.3  Numerical  Simulation  and  Comparison  with  Test  A 

T3q)ical  experimental  results  from  Case  1  of  Test  A  (see  Fig.  31)  are  replotted 
in  Fig.  36,  with  RMS  force  and  displacements  as  a  function  of  wheel  surface 
velocity  (m/s).  Several  speed  regions  were  found  by  measurement  and 
observation.  (1)  When  the  wheel  speed  was  lower  than  6  m/s,  the  vehicle  was 
unstable,  with  high-frequency  (flutter)  due  to  a  large  drag  force.  (2)  When  the 
speed  was  between  6  and  13  m/s,  the  vehicle  was  stable  and  amplitudes  of 
displacement  and  force  were  relatively  small.  (3)  When  the  speed  was  between  13 
and  30  m/s,  the  vehicle  was  unstable  and  dominated  by  slide  and  yaw  instability. 
(4)  When  the  speed  was  over  30  m/s,  the  vehicle  was  unstable  in  almost  every 
direction  and  the  oscillation  amplitudes  in  five  directions  increased  significantly. 
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Velocity,  m/s 

Fig.  36.  Experimental  results  from  Case  1  in  Test  A 

Numerical  simulation  of  the  dynamic  response  of  a  five-DOF  maglev  vehicle 
was  carried  out  by  using  vehicle  dimension  input  as  described  in  Fig.  30  and 
Section  3.  Figure  37  shows  the  simulation  results  of  lateral  and  vertical 
displacements  of  vehicle  mass  center  for  various  wheel  speeds,  with  initial 
perturbations  of  zq  and  yo  equal  to  0.1  mm.  The  system  damping  ratio  in  the 
simulation  was  assumed  to  be  2.5%.  For  vertical  displacement,  the  positive  value 
means  the  vehicle  is  moving  toward  the  guideway.  The  following  conclusions  can 
be  drawn.  (1)  When  vehicle  velocity  is  15  m/s,  the  vehicle  is  stable.  Its  amplitudes 
of  vertical  and  lateral  motion  decrease  as  time  increases.  (2)  For  vehicle  velocities 
of  20,  25,  and  30  m/s,  both  lateral  and  vertical  oscillation  amplitudes  increase 
until  the  vehicle  hits  the  guideway.  (3)  With  defined  vehicle  parameters,  vehicle 
motion  is  dominated  by  slide  and  yaw  instabilities,  even  though  heave  instability 
is  still  present.  This  is  caused  by  weak  guidance  force.  (4)  Oscillation  frequency 
varies  with  the  vehicle  velocity  because  of  lift  and  drag  forces.  In  stable  regions, 
frequency  is  quite  low.  (5)  The  current  computer  program  is  unable  to  deal  with 
nonlinearity  when  the  vehicle  hits  the  guideway.  For  all  calculations  of  time 
history,  the  program  automatically  stops  once  any  side  of  a  vehicle  hits  the 
guideway;  therefore,  it  cannot  be  used  to  calculate  RMS  values  of  displacement 
and  power  spectral  density  (PSD)  for  comparison  with  the  experimental  data 
shown  in  Fig.  36. 

Figure  38  shows  time  histories  of  both  lateral  and  vertical  motion  of  a  vehicle 
with  a  velocity  of  30.0  m/s.  Lateral  motion  is  rapidly  developed  and  finally  the 
vehicle  hits  the  guideway.  Notice  that  the  oscillation  period  changes  when  the 
amplitude  of  lateral  motion  increases.  After  the  vehicle  contacts  the  side  wall  of 
an  L-shaped  guideway  (t  >  1.0  s),  the  lateral  motion  of  vehicle  mass  center 
increases  significantly  and  period  is  reduced. 
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Fig.  37.  Simulation  results  of  (a)  lateral 
and  (b)  mass  center  displacements 
in  Test  A  at  various  vehicle  speeds 


Figure  39  shows  five  vehicle  motions,  slip,  heave,  yaw,  pitch,  and  roll,  when 
vehicle  velocity  =  20.0  m/s  and  initial  perturbation  is  0.1  mm  in  slip  and  heave 
motion  and  0.1°  in  pitch  and  roll  motions.  In  this  particular  case,  pitch  and  roll 
motion  is  stable,  while  slip,  heave,  and  yaw  motion  is  unstable.  The  frequency 
and  phase  of  the  five  motions  are  different,  which  may  be  attributed  to  weak 
coupling  of  motion  in  different  directions.  The  slip  and  yaw  motions  are  strongly 
coupled,  but  only  weakly  coupled  to  the  vertical  modes. 
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Fig.  38.  Simulated  time  histories  of  lateral 
and  vertical  displacement  of 
vehicle  that  contacts  side  wall  of 
guideway 


Fig.  39.  Simulation  results  of  displace¬ 
ments  and  rotations  in  Test  A 
with  vehicle  speed  =  20  mis 


4.4  Numerical  Simulation  and  Comparison  with  Test  B 

Typical  results  from  Test  B  (see  Fig.  34)  are  replotted  in  Fig.  40,  with  RMS 
values  of  heave,  slide,  yaw,  pitch,  and  roll  motion  of  the  vehicle  as  a  function  of 
wheel  surface  velocity  ranging  from  15  to  31  m/s,  with  the  lateral  gap  between 
guidance  magnets  and  the  vertical  part  of  the  guide  way  L  equal  to  13.5  mm. 
Large  motion  occurs  between  19  and  23  m/s.  When  the  velocity  is  lower  than 
18  m/s  or  higher  than  24  m/s,  the  vehicle  was  very  stable  (Fig.  40). 
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Velocity,  m/s 

Fig.  40.  Experimental  results  with 
guidance  magnet  gap  =  13.5  mm 
in  Test  B 

Numerical  simulation  was  carried  out  for  the  vehicle  in  Test  B.  The  vehicle 
configuration  and  magnet's  position  were  measured  and  input  into  the  program. 
Detailed  vehicle  dimensions  and  other  parameters  have  been  described  in  Section 
3.2.  In  the  simulation,  the  initial  gaps  between  magnets  and  the  guideway  L  were 
set  as  follows:  lateral  gap,  13.5  mm;  vertical  gap,  equal  to  the  value  at  which  the 
lift  forces  of  four  lift  magnets  were  sufficient  to  balance  the  vehicle  weight  at  given 
wheel  speeds;  initial  perturbations  zq  and  yo  at  vertical  and  lateral  directions, 
0.1  mm;  and  system  damping  ratio,  2.5%.  Note  that  for  vertical  displacement,  the 
positive  value  means  the  vehicle  is  moving  toward  the  guideway. 

Figures  41-44  show  the  time  histories  of  simulation  results  with  vehicle 
velocities  of  15,  20,  25,  and  30  m/s,  respectively.  The  following  conclusions  can  be 
drawn.  (1)  With  strong  guidance  forces  provided  by  four  guidance  magnets,  slide 
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Fig.  41.  Simulation  results  of  displace¬ 
ment  and  rotation  in  Test  B  with 
vehicle  speed  -  15  mis 


motion  decays  very  quickly  (<4  s,  which  seems  to  agree  well  with  experimental 
observations),  and  no  slide  or  yaw  instabilities  are  present  in  whole  speed  range. 
When  the  vehicle  velocity  increases,  the  period  for  slide  motion  to  die  out 
decreases  because  guidance  forces  increase  with  wheel  velocity.  This  conclusion 
can  also  be  applied  to  all  of  the  rotation  motions.  (2)  Pitch  oscillations  are  quite 
small  because,  in  the  calculation,  drag  force  moments  are  balanced  mostly  by  the 
moment  of  the  tether.  However,  because  of  the  rest  of  the  unbalanced  drag 
moment,  some  pitch  angle  offset  is  still  present  when  oscillation  disappears. 
(3)  Oscillating  frequencies  of  various  motions  are  quite  different  for  this 
configuration.  (4)  The  vehicle  motion  is  dominated  by  heave  motion,  and  heave 
instability  occurs  when  vehicle  velocity  is  20  m/s,  which  is  in  a  good  agreement 
with  the  experimental  result  shown  in  Fig.  40. 


43 


Time,  s 

Fig.  42.  Simulation  results  of  displace¬ 
ment  and  rotation  in  Test  B  with 
vehicle  speed  =  20  ml s 

The  drag  force  produced  by  four  lift  magnets  may  play  an  important  role  in 
these  experiments.  Figure  45  shows  the  guidance  force  (actually,  lift  force  to  this 
magnet,  see  Section  2.3,  but  guidance  force  to  the  vehicle)  and  drag  force  of  the 
guidance  magnet  as  a  function  of  vehicle  speed.  Drag  force  is  small  and  remains 
almost  constant  in  the  test  speed  region.  Thus,  the  influence  of  the  drag  force  of 
the  guidance  magnet  on  vehicle  motion  could  be  ignored.  Figure  46  gives  the 
initial  vertical  gap  at  which  magnets  support  the  vehicle,  and  the  drag  force  that 
is  produced  by  the  lift  magnet  at  various  speeds.  As  the  velocity  increases,  the 
vertical  gap  between  the  magnet  and  sheet  guideway  increases  because  the  lift 
force  needed  to  balance  the  vehicle  weight  remains  the  same;  while  in  this  same 
situation,  drag  force  decreases.  This  property  may  affect  not  only  slide,  yaw,  roll, 
and  pitch  motions  as  their  time  to  die  out  is  reduced  with  increased  velocity,  but 
also  heave  motion,  as  seen  in  Figs.  41-44. 


Magnetic  forces  are  critically  important  to  the  analysis  of  any  vehicle 
dynamics,  design  of  guideway  structures  and  fastening,  and  prediction  of  ride 
quality.  These  force  components  are  considered  from  the  standpoint  of  vehicle 
stability.  This  report  documents  magnetic-force  data  obtained  from  both 
experiments  and  calculations  for  a  permanent  magnet  moving  over  a  sheet 
guideway.  All  of  the  experimental  data  have  been  incorporated  into  an  existing 
computer  program  for  simulation  of  nonlinear  dynamics  of  maglev  vehicles.  The 
magnetic  force  data,  which  have  never  been  published,  will  provide  potential 
references  for  future  maglev  programs. 
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Fig.  44.  Simulation  results  of  displace¬ 
ment  and  rotation  in  Test  B  with 
vehicle  speed  =  30  mis 


Fig.  45.  Comparison  of  vehicle  guidance 
and  drag  forces  as  a  function  of 
velocity 
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Fig.  46.  Comparison  of  vertical  gap  and 
drag  force  as  a  function  of  velocity 
when  vehicle  is  at  equilibrium 


Motion-dependent  magnetic  forces  are  the  key  elements  in  modeling  and 
understanding  dynamic  instabilities  of  maglev  systems.  At  this  time,  it  appears 
that  very  limited  data  are  available  for  motion-dependent  magnetic  forces.  Efforts 
will  be  made  to  compile  analytical  results  and  experimental  data  for  motion- 
dependent  magnetic  forces.  When  this  work  is  completed,  recommendations  will 
be  presented  on  the  needs  of  research  on  magnetic  forces.  In  addition,  specific 
methods  to  obtain  motion-dependent  magnetic  forces  will  be  described  in  detail. 

Two  series  of  extensive  experimental  investigations  on  dynamic  stability  of 
maglev  systems  were  conducted  with  a  free  vehicle  moving  on  a  double  L-shaped 
aluminum  guideway  moimted  on  the  top  of  a  rotating  wheel.  Five  modes  (vertical 
heave,  lateral  slip,  pitch,  yaw,  and  roll)  of  the  vehicle  motion  were  measured  in 
experiments  when  the  rotating  speed  of  the  wheel  was  varied.  Instabilities  of  an 
EDS-type  maglev  system  have  been  observed  through  the  experiments.  Stable  and 
unstable  motion  of  the  maglev  vehicle  was  observed  and  recorded.  An  integrated 
experimental/analytical  study  of  stability  characteristics  is  definitely  an 
important  aspect  of  maglev  research  and  must  be  considered  in  the  development 
of  all  maglev  systems. 

Various  methods  can  be  used  to  stabilize  a  maglev  system:  passive 
electrodynamic  primary  suspension  damping,  active  electrodynamic  primary 
suspension  damping,  passive  mechanical  secondary  suspension,  and  active 
mechanical  secondary  suspension.  With  a  better  understanding  of  vehicle 
stability  characteristics,  a  better  control  law  can  be  adopted  to  ensure  a  high  level 
of  ride  comfort  and  safety. 
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Computer  programs  are  needed  to  screen  new  system  concepts,  evaluate 
various  designs,  and  predict  vehicle  response.  A  computer  code  for  nonlinear 
dynamic  simulation  of  maglev  systems  was  extensively  developed  in  this  study.  It 
can  simulate  nonlinear  dynamic  response  of  maglev  systems  with  six  DOF  when 
the  user  inputs  vehicle  and  track  configurations.  It  provides  more  functionality 
for  broad  vehicle/guideway  designs  than  any  previous  codes.  Numerical 
simulations  presented  in  this  study  for  a  five-DOF  maglev  system  are  in  good 
agreement  with  experimental  results.  However,  it  appears  that  the  stability 
characteristics  of  maglev  vehicles  under  different  conditions  have  not  been 
studied  in  detail.  When  information  on  motion-dependent  magnetic  forces 
becomes  available,  the  computer  codes  can  be  improved  significantly. 

The  experimental  and  theoretical  analyses  developed  in  this  study  identify 
basic  stability  characteristics  and  future  research  needs  of  maglev  systems. 
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Appendix _ 

Measured  quasistatic  magnetic- force  data  for  the  magnet  (50.8  x  25.4  x 
6.35  mm),  corresponding  to  Figs.  4-9,  and  the  magnet  (50.8  x  12.7  x  6.35  mm), 
corresponding  to  Fig.  28,  are  listed  in  Tables  A.1-A.13. 
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Table  A.I.  Measured  lift  magnetic  force  (N)  when  Y  =  5.0  mm 


Height,  V=4.55  V  =  7.89  V=  10.62  V=  13.32  V=  18.20  V  =  24.15  V  =  30.10  V  =  36.10 
m _ m/ s _ m/ s _ m/s _ m/ s _ m/ s _ m/s _ m/s  m/s 


5.0000 

1.0451 

2.5957 

6.0000 

0.91766 

2.2373 

7.0000 

0.79077 

1.9308 

8.0000 

0.69039 

1.6812 

9.0000 

0.61040 

1.4675 

10.000 

0.53812 

1.2908 

12.700 

0.38141 

0.91256 

13.000 

0.37110 

0.88029 

15.000 

20.000 

25.000 

30.000 

35.000 

40.000 

45.000 

50.000 

0.28903 

0.69668 

3.9378 

5.2324 

7.2951 

3.3923 

4.5002 

6.2764 

2.9271 

3.8891 

5.4203 

2.5420 

3.3842 

4.7035 

2.2207 

2.9615 

4.1118 

1.9460 

2.5938 

3.6103 

1.3865 

1.8570 

2.5750 

1.3383 

1.7943 

2.4931 

1.0637 

1.4284 

1.9862 

9.2438 

10.562 

11.651 

7.9298 

9.0468 

9.9822 

6.8365 

7.7943 

8.5869 

5.9278 

6.7518 

7.4294 

5.1641 

5.8795 

6.4625 

4.5277 

5.1504 

5.6501 

3.2354 

3.6788 

4.0371 

3.1303 

3.5554 

3.9008 

2.4957 

2.8386 

3.1147 

1.8988 

1.2636 

0.90294 

0.69926 

0.57776 

0.48013 

0.36727 
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Table  A.2.  Measured  lift  magnetic  force  (N)  when  Y  -  12.7  mm 


Height, 

mm 

5.0000 

6.0000 

7.0000 

8.0000 

9.0000 

10.000 

12.700 

13.000 

15.000 

20.000 

25.000 

30.000 

35.000 

40.000 

45.000 

50.000 


V  =  4.55 
m/ s 
1.0955 
0.94664 
0.81583 
0.70712 
0.61203 
0.53374 
0.36982 
0.35219 
0.27539 


V  =  7.89 
m/s 
2.6161 
2.2545 
1.9367 
1.6686 
1.4513 
1.2658 
0.88409 
0.84933 
0.66572 


V  =  10.62 
m/s 
3.9605 
3.3970 
2.9226 
2.5221 
2.1937 
1.9114 
1.3357 
1.2875 
1.0054 


V=  13.32 
m/s 
5.2814 
4.5139 
3.8977 
3.3668 
2.9296 
2.5579 
1.7984 
1.7361 
1.3636 


V=  18.20 
m/ s 
7.3515 
6.3151 
5.4191 
4.6907 
4.0759 
3.5538 
2.4904 
2.4037 
1.8837 


V  =  24.15 
m/s 
9.3087 
7.9683 
6.8397 
5.8954 
5.1076 
4.4472 
3.1052 
2.9957 
2.3421 


V  =  30.10 
m/ s 
10.676 
9.1351 
7.8049 
6.7165 
5.8145 
5.0512 
3.5267 
3.4014 
2.6608 


V  =  36.10 
m/s 
11.805 
10.080 
8.6117 
7.3991 
6.3945 
5.5470 
3.8584 
3.7259 
2.9041 
1.6615 
1.0252 
0.66869 
0.45884 
0.32739 
0.23401 
0.16809 
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Table  A.3.  Measured  lift  magnetic  force  (N)  when  h  -  7.0  mm 


Gap, 

mm 

V=4.55 

m/s 

V=7.89 

m/s 

V=  10.62 
m/s 

V=  13.32 

m/s 

V=  18.20 

m/s 

V=24.15 

m/s 

V  =  30.10 
m/s 

V=  36.10 

m/s 

5.0000 

0.80121 

1.9141 

2.9061 

3.8765 

5.4086 

6.8017 

7.7810 

8.5608 

6.0000 

0.80485 

1.9186 

2.9103 

3.8865 

5.4090 

6.7890 

7.7679 

8.5471 

7.0000 

0.80143 

1.9207 

2.9116 

3.8811 

5.4047 

6.7887 

7.7602 

8.5383 

8.0000 

0.80612 

1.9224 

2.9143 

3.8877 

5.4072 

6.7844 

7.7642 

8.5411 

9.0000 

0.80847 

1.9266 

2.9177 

3.8896 

5.4110 

6.7777 

7.7583 

8.5448 

10.000 

0.81084 

1.9339 

2.9212 

3.8907 

5.4181 

6.7858 

7.7670 

8.5598 

12.700 

0.80694 

1.9314 

2.9187 

3.8878 

5.4097 

6.8048 

7.7962 

8.5949 

13.000 

0.80654 

1.9298 

2.9218 

3.8918 

5.4166 

6.8056 

7.7969 

8.6008 

15.000 

17.000 

0.80960 

1.9306 

2.9161 

3.8883 

5.4222 

6.8127 

7.8057 

7.8208 

8.6321 

8.6663 

20.000 

22.000 

25.000 

26.000 

27.000 

28.000 

29.000 

30.000 

31.000 

32.000 

35.000 

40.000 

45.000 

50.000 


8.7054 

8.7193 

8.7220 

8.7130 

8.7019 

8.6783 

8.6503 

8.6122 

8.5553 

8.4845 

8.1381 

6.9166 

4.9871 

2.7338 
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Table  A.4.  Measured  lift  magnetic  force  (N)  when  h  -  12.7  mm 


Gap, 

m  m 

V=4.55 
m/ s 

V=7.89 

m/s 

5.0000 

0.37788 

0.90238 

6.0000 

0.37967 

0.90227 

7.0000 

0.37603 

0.90000 

8.0000 

0.37182 

0.89802 

9.0000 

0.37459 

0.89189 

10.000 

0.37172 

0.88857 

12.700 

0.36636 

0.88496 

13.000 

0.36895 

0.88526 

15.000 

0.36174 

0.87754 

20.000 

25.000 

30.000 

35.000 

40.000 

45.000 

50.000 


V=  10.62 

m/s 

V  =  13.32 

m/s 

V=  18.20 

m/s 

1.3718 

1.8505 

2.5732 

1.3676 

1.8437 

2.5494 

1.3620 

1.8343 

2.5468 

1.3575 

1.8266 

2.5321 

1.3531 

1.8208 

2.5151 

1.3477 

1.8163 

2.5118 

1.3382 

1.7973 

2.4896 

1.3369 

1.7979 

2.4892 

1.3247 

1.7828 

2.4769 

V  =  24.15 

m/ s 

V  =  30.10 

m/s 

V  =  36.10 

m/s 

3.2065 

3.6706 

4.0127 

3.1774 

3.6340 

3.9705 

3.1599 

3.6043 

3.9352 

3.1395 

3.5848 

3.9074 

3.1268 

3.5595 

3.8855 

3.1097 

3.5498 

3.8729 

3.0842 

3.5213 

3.8490 

3.0826 

3.5163 

3.8488 

3.0673 

3.5092 

3.8402 

3.8250 

3.7733 

3.6226 

3.2929 

2.7036 

1.8991 

1.0999 
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Table  A.  5.  Measured  guidance  magnetic  force  (N)  when  Y  =  5.0  mm 


Height,  V=4.55  V=7.89  V=  10.62  V=  13.32  V=  18.20  V  =  24.15  V  =  30.10  V=  36.10 

mm  m/s  m/s  m/s  m/s _ m/s _ m/s _ m/ s _ m/s 


5.0000  0.18065  0.44393 

6.0000  0.18417  0.45186 

7.0000  0.18068  0.46084 

8.0000  0.18766  0.46147 

9.0000  0.18833  0.47170 

10.000  0.18729  0.47215 

13.000  0.19191  0.47090 

15.000  0.18777  0.47444 

20.000 
25.000 
30.000 
35.000 
40.000 


0.71410 

0.98960 

1.5682 

0.72744 

1.0170 

1.6034 

0.74301 

1.0289 

1.6325 

0.74814 

1.0445 

1.6458 

0.75691 

1.0496 

1.6620 

0.75790 

1.0510 

1.6666 

0.76569 

1.0615 

1.6699 

0.76355 

1.0576 

1.6677 

2.2220 

2.7983 

3.3500 

2.2619 

2.8586 

3.4133 

2.2890 

2.8935 

3.4652 

2.3130 

2.9141 

3.4866 

2.3285 

2.9345 

3.5110 

2.3385 

2.9420 

3.5156 

2.3323 

2.9421 

3.5116 

2.3194 

2.9147 

3.4810 

3.3748 

3.2648 

3.1568 

3.0531 

2.9309 
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Table  A.  6.  Measured  guidance  magnetic  force  (N)  when  Y  =  12.7  mm 


Height, 

m  m 

V=4.55 

m/ s 

V=7.89 

m/s 

V=  10.62 

m/s 

V=  13.32 

m/s 

V=  18.20 

m/s 

V  =  24.15 

m/s 

V  =  30.10 

m/s 

V=  36.10 

m/s 

5.0000 

0.097706 

0.20690 

0.30712 

0.38804 

0.53884 

0.65595 

0.74890 

0.82918 

6.0000 

0.094964 

0.20529 

0.30992 

0.39729 

0.56286 

0.69580 

0.80601 

0.90144 

7.0000 

0.091036 

0.20977 

0.31310 

0.40684 

0.57802 

0.73006 

0.84524 

0.95586 

8.0000 

0.090508 

0.20402 

0.31586 

0.41019 

0.59173 

0.75649 

0.88238 

0.99907 

9.0000 

0.093734 

0.20588 

0.31544 

0.41667 

0.60427 

0.77338 

0.91097 

1.0331 

10.000 

0.087470 

0.20646 

0.31930 

0.42212 

0.61046 

0.78828 

0.93361 

1.0681 

12.700 

0.083869 

0.19930 

0.30894 

0.41791 

0.61670 

0.80983 

0.96883 

1.1137 

13.000 

0.081553 

0.19838 

0.31017 

0.41911 

0.61545 

0.81066 

0.97017 

1.1160 

15.000 

0.079007 

0.19274 

0.30440 

0.41200 

0.61152 

0.81318 

0.97188 

1.1252 

20.000 

25.000 

30.000 

35.000 

40.000 

45.000 

50.000 


1.1000 

1.0424 

0.96882 

0.89044 

0.82608 

0.76790 

0.72454 
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Table  A.7.  Measured  guidance  magnetic  force  (N)  when  h  =  7.0  mm 


Gap, 

m  m 

V  =  4.55 

m/s 

V  =  7.89 

m/s 

V=  10.62 
m/s 

5.0000 

0.19122 

0.45731 

0.73137 

6.0000 

0.16639 

0.39968 

0.63375 

7.0000 

0.15097 

0.34776 

0.55057 

8.0000 

0.13917 

0.31218 

0.48647 

9.0000 

0.13029 

0.28606 

0.43460 

10.000 

0.12170 

0.26143 

0.39458 

13.000 

0.11148 

0.22066 

0.32044 

15.000 

0.10777 

0.21193 

0.30006 

16.000 

17.000 

18.000 

20.000 

22.000 

25.000 

23.000 

22.000 

24.000 

25.000 

26.000 

27.000 

28.000 

29.000 

30.000 

31.000 

32.000 

33.000 

34.000 

35.000 

36.000 

37.000 

38.000 

39.000 

40.000 

41.000 

42.000 

43.000 


V=  13.32  V=  18.20  V=24.15  V  =  30.10  V=  36.10 


m/s 

m/s 

m/s 

m/s 

m/s 

1.0322 

1.6304 

2.2933 

2.8992 

3.4594 

0.88455 

1.3820 

1.9332 

2.4178 

2.8615 

0.76025 

1.1773 

1.6356 

2.0328 

2.3895 

0.66585 

1.0152 

1.3904 

1.7220 

2.0064 

0.58761 

0.88464 

1.1986 

1.4675 

1.6938 

0.52755 

0.77599 

1.0376 

1.2542 

1.4431 

0.41785 

0.56402 

0.71816 

0.83196 

0.92319 

0.37833 

0.48031 

0.59664 

0.66109 

0.71464 

0.62420 

0.56259 

0.51805 

0.47121 

0.45942 

0.53118 

0.47426 

0.46136 

0.50279 

0.54166 

0.57850 

0.63566 

0.70825 

0.78758 

0.87856 

0.97932 

1.1058 

1.2436 

1.3905 

1.5497 

1.7205 

1.8941 

2.0679 

2.2402 

2.3956 

2.5338 

2.6490 

2.7376 


eo 


Table  A.  7.  Cont  'd 


Gap,  V=4.55  V  =  7.89  V=  10.62  V=  13.32  V=  18.20  V  =  24.15  V  =  30.10  V=:  36.10 
mm  m/s  m/s _ m/s _ m/s _ m/s _ m/ s _ va!  s _ 


44.000 

2.7861 

45.000 

2.8098 

46.000 

2.7861 

46.000 

2.7877 

47.000 

2.7349 

48.000 

2.6488 

49.000 

2.5441 

50.000 

2.4179 

51.000 

2.2653 

52.000 

2.1055 

53.000 

1.9357 
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Table  A.8.  Measured  guidance  magnetic  force  (N)  when  h  =  12.7  mm 


Gap, 

m  m 

V=4.55 

m/s 

V  =  7.89 

m/s 

V=  10.62 

m/s 

V=  13.32 
m/s 

V=  18.20 

m/s 

V  =  24.15 

m/s 

V  =  30.10 
m/s 

V  =  36.10 

m/s 

0.16979 

0.45445 

0.74629 

1.0482 

1.6583 

2.3402 

2.9856 

3.5008 

0.15160 

0.38587 

0.64142 

0.89846 

1.4110 

1.9835 

2.5116 

2.9385 

0.13681 

0.33966 

0.56009 

0.77435 

1.2130 

1.6944 

2.1361 

2.4894 

MiU 

0.11769 

0.30203 

0.49382 

0.67819 

1.0556 

1.4616 

1.8246 

2.1221 

9.0000 

0.10520 

0.26779 

0.43972 

0.60001 

0.92745 

1.2662 

1.5761 

1.8277 

10.000 

0.094461 

0.24160 

0.38845 

0.53316 

0.81830 

1.1116 

1.3697 

1.5795 

12.700 

0.078118 

0.19241 

0.31376 

0.41773 

0.61378 

0.81021 

0.97368 

1.1036 

13.000 

0.078038 

0.18856 

0.30458 

0.40892 

0.59571 

0.78456 

0.94053 

1.0712 

15.000 

0.074923 

0.17032 

0.27774 

0.36331 

0.51582 

0.65495 

0.76163 

0.85348 

20.000 

25.000 

26.000 

27.000 

28.000 

29.000 

30.000 

31.000 

32.000 

33.000 

34.000 

35.000 

38.000 

40.000 

42.000 

44.000 

45.000 

46.000 

47.000 

48.000 

49.000 

50.000 


0.59110 

0.56213 

0.58183 

0.60261 

0.62805 

0.65752 

0.69215 

0.73462 

0.77862 

0.82802 

0.87711 

0.92915 

1.0834 

1.1639 

1.2147 

1.2380 

1.2343 

1.2158 

1.1890 

1.1532 

1.1144 

1.0627 
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Table  A.  9.  Measured  drag  magnetic  force  (N)  when  Y  -  5.0  mm 


Height, 
m  m 


V=4.55 
m/ s 


y  _  7  39  V  -  10.62  V  =  13.32  V  =  18.20  V  =  24.15  V  =  30.10  V  =  36.10 

m/s  m/ s _ m/s _ m/s _ m/ s _ m/s _ 


5.0000 

3.0899 

4.5607 

6.0000 

2.7156 

4.0115 

7.0000 

2.4139 

3.5401 

8.0000 

2.1504 

3.1631 

9.0000 

1.9392 

2.8430 

10.000 

1.7548 

2.5734 

12.700 

1.3775 

2.0225 

13.000 

1.3439 

1.9703 

15.000 

1.1535 

1.7049 

20.000 

25.000 

30.000 

35.000 

40.000 

45.000 

50.000 

5.2778 

5.6976 

5.9508 

4.6392 

5.0186 

5.2514 

4.1140 

4.4653 

4.6733 

3.6723 

3.9933 

4.1989 

3.3110 

3.6114 

3.8042 

3.0068 

3.2859 

3.4736 

2.3846 

2.6338 

2.7998 

2.3329 

2.5755 

2.7420 

2.0209 

2.2399 

2.4076 

5.8713 

5.6439 

5.4027 

5.1931 

5.0020 

4.7839 

4.6290 

4.4703 

4.2780 

4.1691 

4.0321 

3.8626 

3.7907 

3.6719 

3.5230 

3.4708 

3.3689 

3.2403 

2.8228 

2.7617 

2.6715 

2.7733 

2.7178 

2.6260 

2.4519 

2.4191 

2.3380 

1.8856 

1.6287 

1.4599 

1.3414 

1.2395 

1.1352 

1.0074 
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Table  A.  10.  Measured  drag  magnetic  force  (N)  when  Y  =  12.7  mm 


Height,  V=4.55  V=7.89  V=  10.62  V=  13.32  V=  18.20  V  =  24.15  V  =  30.10  V=  36.10 

mm _ m/s _ m/s _ m/s  m/s  m/s  m/s  m/s  m/s 


5.0000 

2.7672 

3.9798 

6.0000 

2.3707 

3.4231 

7.0000 

2.0555 

2.9609 

8.0000 

1.7866 

2.5791 

9.0000 

1.5689 

2.2634 

10.000 

1.3802 

1.9988 

12.700 

1.0130 

1.4573 

13.000 

0.98104 

1.4158 

15.000 

0.79730 

1.1534 

20.000 

25.000 

30.000 

35.000 

40.000 

45.000 

50.000 

4.5695 

4.9081 

5.0506 

3.9263 

4.2150 

4.3393 

3.3931 

3.6459 

3.7496 

2.9554 

3.1830 

3.2622 

2.5962 

2.7984 

2.8714 

2.2902 

2.4781 

2.5396 

1.6869 

1.8307 

1.8816 

1.6327 

1.7756 

1.8265 

1.3391 

1.4644 

1.5105 

4.9185 

4.7044 

4.4699 

4.2126 

4.0277 

3.8128 

3.6363 

3.4675 

3.2863 

3.1675 

3.0208 

2.8608 

2.7839 

2.6552 

2.5075 

2.4621 

2.3552 

2.2176 

1.8267 

1.7532 

1.6606 

1.7807 

1.7090 

1.6166 

1.4749 

1.4251 

1.3461 

0.93444 

0.71183 

0.57532 

0.47564 

0.40846 

0.34486 

0.29709 
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Table  A.  12.  Measured  drag  magnetic  force  (N)  when  h  =  12.7  mm 


Gap,  V=4.55  V=7.89  V=  10.62  V=  13.32  V=  18.20  V  =  24.15  V  =  30.10  V  =  36.10 

mm _ m/s _ m/s _ m/s  m/s  m/s  m/s  m/s  m/s 


1.3715 

2.0248 

1.2988 

1.9054 

1.2261 

1.7999 

8.0000 

1.1730 

1.7116 

9.0000 

1.1283 

1.6430 

10.000 

1.0857 

1.5840 

12.700 

1.0074 

1.4657 

13.000 

1.0024 

1.4542 

15.000 

0.96331 

1.3938 

20.000 

25.000 

30.000 

35.000 

40.000 

45.000 

50.000 

2.3827 

2.6163 

2.7915 

2.2245 

2.4369 

2.5863 

2.0951 

2.2932 

2.4162 

1.9946 

2.1750 

2.2757 

1.9075 

2.0744 

2.1626 

1.8357 

1.9908 

2.0695 

1.6939 

1.8303 

1.8864 

1.6756 

1.8171 

1.8685 

1.6080 

1.7374 

1.7872 

2.8030 

2.7596 

2.6502 

2.5708 

2.5234 

2.4080 

2.3824 

2.3231 

2.2141 

2.2347 

2.1749 

2.0624 

2.1101 

2.0519 

1.9459 

2.0165 

1.9469 

1.8411 

1.8246 

1.7553 

1.6566 

1.8045 

1.7324 

1.6420 

1.7157 

1.6522 

1.5528 

1.4401 

1.3995 

1.3726 

1.3360 

1.2283 

1.0341 

0.77217 


Table  A.  13.  Measured  lift  and  drag  magnetic  force  (N)  for  guidance  magnet 


Velocity, 

m/s 

Lift,  h  = 
8.5  mm 

Lift,  h  = 
13.5  mm 

Lift,  h  = 
18.5  mm 

Drag,  h  = 
8.5  mm 

Drag,  h  = 
13.5  mm 

Drag,  h  = 
18,5  mm 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

2.0000 

0.030830 

0.015949 

0.0090586 

0.18331 

0.091968 

0.051258 

4.0000 

0.11647 

0.059970 

0.033941 

0.34709 

0.17321 

0.096164 

6.0000 

0.24052 

0.12299 

0.069256 

0.47963 

0.23748 

0.13110 

8.0000 

0.38513 

0.19534 

0.10934 

0.57870 

0.28387 

0.15566 

10.000 

0.53638 

0.26971 

0.15002 

0.64827 

0.31484 

0.17142 

12.000 

0.68549 

0.34172 

0.18889 

0.69452 

0.33391 

0.18051 

14.000 

0.82778 

0.40921 

0.22484 

0.72344 

0.34439 

0.18489 

16.000 

0.96118 

0.47138 

0.25753 

0.73993 

0.34890 

0.18607 

18.000 

1.0852 

0.52820 

0.28704 

0.74769 

0.34938 

0.18516 

20.000 

1.1999 

0.57995 

0.31361 

0.74940 

0.34718 

0.18290 

22.000 

1.3061 

0.62709 

0.33755 

0.74696 

0.34325 

0.17982 

24.000 

1.4043 

0.67010 

0.35916 

0.74172 

0.33823 

0.17626 

26.000 

1.4953 

0.70943 

0.37873 

0.73462 

0.33256 

0.17245 

28.000 

1.5799 

0.74551 

0.39652 

0.72633 

0.32656 

0.16854 

30.000 

1.6587 

0.77871 

0.41274 

0.71731 

0.32041 

0.16464 

32.000 

1.7323 

0.80935 

0.42760 

0.70789 

0.31426 

0.16080 

34.000 

1.8011 

0.83773 

0.44124 

0.69831 

0.30820 

0.15708 

36.000 

1.8657 

0.86408 

0.45383 

0.68872 

0.30229 

0.15348 

38.000 

1.9265 

0.88863 

0.46547 

0.67923 

0.29656 

0.15004 

40.000 

1.9837 

0.91156 

0.47628 

0.66993 

0.29104 

0.14675 
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